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0^ ■ We show that the physical mechanism of population transfer in a 3-level system with a closed 

loop of coherent couplings (loop-STIRAP) is not equivalent to an adiabatic rotation of the dark-state 
^ , of the Hamiltonian but coresponds to a rotation of a higher-order trapping state in a generalized 

adiabatic basis. The concept of generalized adiabatic basis sets is used as a constructive tool to 
design pulse sequences for stimulated Raman adiabatic passage (STIRAP) which give maximum 
population transfer also under conditions when the usual condition of adiabaticty is only poorly 
fulfillcd. Under certain conditions for the pulses (generalized matched pulses) there exists a higher- 
order trapping state, which is an exact constant of motion and analytic solutions for the atòmic 
dynamics can be derived. 

o ■ 
o ■ 

qq , I. INTRODUCTION 
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The transfer of population in a multi-level atòmic system from an initial to a target quantum state in a fast and 
f ï effective way is currently a problem of practical importance as well as of substantial theoretical interest. If there is 
^Jj a dipole allowed transition between an initial and a target state, one can achieve the desired transfer by using eithcr 
a constant-frequeney 7r-pulse tuned to resonance, or an adiabatic process based on a swept carrier-frequeney. Since 
a dipole-allowed transition implies radiative decay, one is however often interested in systems with two metastable 
states without a direct electric-dipole coupling. Whereas an extension of the two-state 7r-pulse approach to multistate 
excitation is possible, these techniques require careful control of the pulse areas. Adiabatic processes do not require 
. ' such precise control, if the time-evolution is slow (meaning, generally, large pulse areas). In a three-state Raman- 
transition system, for example, it is possible to achieve adiabatic passage with the use of two constant-frequeney pulses 
suitably delayed (counterintuitive order) The process of this stimulated Raman adiabatic passage (STIRAP) 
can be represented by a slow rotation of a decoupled eigenstate of the H amiltonian (dark state) Q . 

The disadvantage of STIRAP is the requirement for large pulse areas: to ensure adiabatic time evolution the effective 
average Rabi-frequency of the pulses must be large compared to the radiative decay rates of the intermediate level(s). 
Non-adiabatic corrections and the associated diabatic losses |5|-fjJ, scale with 1/flT where 7ïf2 is a characteristic 
interaction energy and T is the effective time required for the transfer. In some potential applications, as for example 
the transfer of information in form of coherences ||, it is desirable to minimize these losses without the need of 
intense pulses or long transfer times. Intense fields induce time-varying ac-Stark shifts, which may be detrimental to 
the coherence transfer. Short times are required to minimize the effect of decoherence processes during the transfer 
§• 

An approach, which reduces non-adiabatic losses for pulses of moderate fluence in a three-state system, was recently 
introduced in Ref. |Ï(J. In addition to the pair of Raman pulses ("pump pulse" and "Stokes pulse") which couplc 
the initial and target state via a common upper level, a direct coupling (called "detuning pulse") between them is 
introduced. This scheme of loop-STIRAP does not require the usual adiabaticity conditions (of large pulse areas), nor 
is it of the 7r-pulse type (requiring specific pulse areas). Nevertheless, the scheme can produce complete population 
transfer. 

In the present paper we show that the physical mechanism of loop-STIRAP is not an adiabatic rotation of the dark 
state, but the rotation of a higher-order trapping state in a generalized adiabatic basis. The concept of generalized 
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adiabatic basis sets allows to rationalize many examples of population transfer even when the adiabaticity condition is 
poorly fulfilled. If pump and Stokes pulses fulfill certain conditions (they are then called generalized matched pulses), 
a higher-order trapping exists, which is an exact constant of motion. In this case analytic solutions for the atòmic 
dynamics can be found which in contrast to the case of ordinary matched pulses with identical pulse shape JO] also 
include the possibility of population transfer. This can be exploited to design pulse sequences which give maximum 
population transfer. In contrast to techniques based on optimum control theory, which are used for such tasks, the 
generalized-dark-state concept provides a physical interpretation of the results. However, the design of pulse, which in 
some cases can lead to complete population transfer (i.e. without any diabatic losses) needs to respect more restrictive 
requirements for specific pulse properties similar to 7r-pulse techniques. 

Our paper is organized as follows. In Sec. II we discuss the loop-STIRAP and propose a simple physical interpretation 
in terms of an adiabatic rotation of a generalized trapping state. In Sec. III we define generalized trapping states via 
an iterative partial diagonalization of the time-dependent Hamiltonian. In Sec. IV we derive conditions under which 
a higher-order trapping state is an exact constant of motion and thus allow for an analytic solution of the atòmic 
dynamics. Finally, various examples of population and coherence transfer based on generalized trapping states are 
discussed in Sec.V. 



II. LOOP-STIRAP 



To set the stage we consider in the present section a three-state system driven by coherent fields in a loop con- 
figuration, as shown in Fig. The bare atòmic states tpi and -03 are coupled by a resonant Raman transition via 
the excited atòmic state 1^2 by a pump pulse and a Stokes pulse, having Rabi- freqüències P{t) and S(t), respectively, 
which are in general complex. In addition there is a direct coupling between states 1 and 3 by a coherent detuning 
pulse described by the (complex) Rabi-frequency D(t). Before the application of the pulses the system is in state 1 
and the goal is to transfer all population into the target state 3 by an appropriate sequence of pulses. For simplicity we 
assume that the carrier freqüències of the pulses coincide with the atòmic transition freqüències and that the phases 
of the pulses are time- independent. Sincc the phases of pump and Stokes fields can be included into the definition of 
the bare atòmic states ipi and ^3, they can be set equal to zero without loss of generality. The phase of the detuning 
pulse is relevant and cannot be eliminated. The time-dependent Schròdinger equation for this system, in the usual 
rotating wave approximation, reads 



-C(t) = -iW(i)C(í) 



(1) 



where C(í) is the column vector of probability amplitudes C n (t) — 
matrix W(í) has the form 



[n\ip(t)), (\n) e {tpi,ip 2 ,ip3})- The evolution 



W(t) = \ 





P(t) 

D*(t) 



P(t) D(t) 
S(t) 
S(t) 



(2) 




FIG. 1. Three-state system with loop linkage. P(t), S(t), D(t) denote Rabi- freqüències of pump, Stokes and detuning pulse. 

It is well known that the counterintuitive pulse sequence (Stokes puls precedes pump pulse, without a detuning 
pulse) leads to an almost complete population transfer, if the adiabaticity condition ÍIT ^> 1 is fulfilled. Here T is 
the characteristic time for the transfer, given by the interval where S(t) and P(t) overlap, and Çl the effective total 
Rabi-frequency averaged over the interval T 
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n = = / dt^/ïW+sW- 

J — oo 



(3) 



As shown in |Ï0[ an almost perfect transfer is also possible when pump and Stokes alone do not fulfill the adiabaticity 
condition by applying an additional detuning pulse. Fig. |^ illustrates an example of ramped pump and Stokes pulses 
intersected by a hyperbolic-secant detuning pulse. 



r 1 7r 

P(t) = A P sin - arctan(í/Tp) + - 

rl 7r" 
S(t) = A s cos - arctan(í/r s ) + - 

D(t) = A D sech[t/T D ]. 



(4) 

(5) 
(6) 
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FIG. 2. Pair of ramped pump (line) and Stokes (dotted) pulses with Ap = As = 20 and Tp = Ts = 0.1 applied in 
counterintuitive order (Stokes precedes pump) with additional hyperbolic secant detuning pulse (dashed) with Ap> = —13 Ai 
and To = 0.2 

Fig. |^ shows examples of population histories for these pulses. When only the pump and Stokes pulses are present, 
the population transfer is rather poor, since the pulse areas are small (Í1T ~ \Ap\ Tp = \As\ Ts — 2). As can be seen 
from the upper part of Fig||, only about 70% of the initial population ends up in state 3. 
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FIG. 3. Populations of states xpi (line), ipa (dotted) and ipz (dashed) for pulse sequence of Fig§ The upper picture shows 
population when only pump and Stokes pulses are applied, and the lower one if the detuning pulse is added. 
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The situation is remarkably different when a detuning pulse with |Ad|Td 2.7 and a phase factor of e - ™/ 2 is 
applied; see the lower part of Fig.||. With a detuning pulse present all the population is transfered from the initial to 
the target state. This result is relatively insensitive to changes in the amplitude (or the shape of the detuning pulse) 
if the phase is — 7r/2. 

We note that in contrast to ordinary STIRAP there is (for a short time) a substantial intermediate population of 
state 2. This indicates that the transfer does not occur as adiabatic rotation of the dark state from ipi to ■03- 

For our present discussion it is useful to describe ordinary STIRAP in tcrms of the following set of adiabatic 
superposition states 



with the unitary matrix 



■$l(t)" 






$a(t) 


= U(t)* 


^2 


_* 3 (t)_ 









1 






(7) 



U(t) = 



sin6» (<) cos 6» (i) 
i cos 6*o (í) — isin6o(t) 



The dynamical angle 8q is defined by 



tan 6o(t) 



S(t)- 



(8) 



(9) 



The vector of probability amplitudes in the bare atòmic basis C(í) and a corresponding vector B(í) in the superposition 
basis (Q) are related through the transformation 



B(í) = U(í)C(í). 

Since U(í) is time-dependent, the transformed evolution matrix has the form 

W(t) -> W(t) - U^W^UW" 1 + i Ú(í)U(t) _1 . 



(10) 



(11) 



In the adiabatic limit, the second term can be disregarded and we are left with the first one, which for ordinary 
STIRAP, i.e. without the detuning pulse, reads 



U(i)W(t)U(í) _1 



Q(t) 
íï(í) 




(12) 



where Q(í) = J P(t) 2 + S(t) 2 . One recognizes that the superposition state $3(í) is decoupled from the coherent inter- 
action in this limit. Moreover, because ^(í) does not contain the excited atòmic state ip2, it does not spontaneously 
radiate and is therefore called a dark state Q. For a counterintuitive sequence of pulses the angle 9o(t) vanishes 
initially and approaches n/2 for t — > oo. Thus 3>3(í) asymptotically coincides with the initial and target states for 
i — > ±oo respectively. Therefore ordinary STIRAP can be understood as a rotation of the adiabatic dark state <I>3 (í) 
from the initial to the target bare atòmic state . Non-adiabatic corrections are contained in the second contribution 
to W(í) 



z Ü(í)U(í)- 1 = 



1 




20 o (í) 
2Ò (t) 



(13) 



They give rise to a coupling between the dark state ^(í) and the so-called bright state $2^)- 

Let us now apply the same transformation to the loop-STIRAP system, i.e. including the detuning pulse. We find: 



W(í) 



o n(t) 

Q(t) Re[D(t)] sin26» (í) 

26 (t) - i [D* (t) sin 2 6 (t)-D (í) cos 2 6 (í)] 





2Ó (t)+i[D(t) sin 2 O (t) -D*(t) cos 2 ( 
-Re[D{t)] sin26» (í) 



)(*)] 



(14) 
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If D(t) is real or complex but not strictly imaginary, there is a time dependent energy shift of the superposition states 
$2(í) and $3(í) and the detuning pulse adds an imaginary part to the nonadiabatic coupling. If D(t) is imaginary, as 
in the example discussed above, there is no detuning but a real^ contribution to the nonadibatic coupling. Let us now 
assume an imaginary detuning pulse, i.e. D(t) — iD(t), with D(t) being real. In this case the transformed evolution 
matrix simplifies to 



W(t) = - 





n(í) 
o 



fi(í) 

o 

2Ò (t) - D(t) 





2Ò (t)-D{t) 




(15) 



If the amplitude of the detuning pulse matches the non-adiabatic term, i.e. if D (i) = 2Òo(t), the dark state <E>3 is 
exactly decoupled even if the adiabaticity condition for pump and Stokes alonc (fi(í) bcing much larger than 0o(t)) 
is not fulfilled. However, since 9o(t) rotates from to tt/2, the detuning pulse would have to be exactly a 7r-pulse in 
such a case. 



dtD(t) = / dt2Ò (t) = 20 o (í) 



+ OO 



(16) 



Furthermore no pump or Stokes pulses were required for population transfer to begin with, since at any time the 
entire population is kept in the dark state by the action of the detuning pulse and thus pump and Stokes would not 
interact with the àtoms. This is consistent with the observation that an exactly decoupled state $3 implies exactly 
vanishing (not only adiabatically small!) probability amplitude of the excited bare state ip2 for all times. Since the 
origin of population transfer in this case is the well-known phenomenon of 7r-pulse coupling, which requires a careful 
control of the area and the shape of the detuning pulse, the case D(t) = 29q{£) is of no further interest here. 

On the other hand, if D(t) is negative, as in the example of Fig.||, the non-adiabatic coupling is effectively increased 
by the detuning pulse (note that d9o(t)/dt > 0). Thus the success of population transfer in Fig.y cannot be understood 
as dark-state rotation. This is illustrated in Fig.^, which shows the populations of the superposition states $1 = -02, 
$2, and $3 for the above example. One clearly sees that about 80% of the population is driven out of the dark state 
during the interaction. 




FIG. 4. Population of superposition states $1 (dashed), $2 (dotted), and the dark state $3 (line). Parameters are that of 
Fig| 

It is worth noting, however, that <ï>2 remains almost unpopulated during the interaction and all population exchange 
happens between states $1 and $3. This suggests an interpretation of the process as adiabatic population return 
between the superposition states $1 and $3. In fact comparing the dressed-state evolution matrix W(í), Eq. (|Ï5|) , with 
the bare-state evolution matrix W(í), Eq.(|^) (without detuning pulse), one recognizes a formal agreement with the 
correspondence P(t) «-> Q(t) and S(t) <-> 29o(t) — D(t). That is there exists a generalized trapping state which is a 
superposition of the states <E>i and $3. Since here Í2(í) = const. and 29o(t) — D(t) vanishes in the asymptotic limits 
t — > ±00, this generalized trapping state coincides with $3 for t — > ±00, which in turn coincides with ipi and ^3 in 
the respective limits. 

To quantify this statement let us introduce a basis of second-order adiabatic states. Using now the first-order states 
$1, $2, and $3 as a basis set instead of the bare atòmic states, we introduce in analogy to Eq.(R) 
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>1 




= Ui(í)* 


$2(í) 


= Ui(t)* ■ U(i)* 


^2 


4 2) (í) 




$3(0 




_^3 



(17) 



The unitary transformation matrix is given by 

Ui(í) = 



1 

sin0i(t) 
i cos 0i (í) 





COS 01 (í) 

-i sin0i(í) 



(18) 



(19) 



which has the same form as U(t), Eq.(g) but here the dynamical angle 0i(í) is defined by 

tan 0i (í) = . . 

20 o (í) - £>(t) 

Denoting the vector of probability amplitudes in these generalized adiabatic states by B^ 2 )(í) we find the relation 

B( 2 )(í) = Ui(í)B(í). (20) 

One easily verifies that for the above example more than 95% of the population remains in the generalized trapping 
state $3 (í). Thus the success of the population transfer in loop STIRAP can be understood as a rotation of the 

(2) 

second-order decoupled state $3 (i) - which is an approximate constant of motion - from the initial to the target 
bare atòmic state. 



III. GENERALIZED ADIABATIC BASIS AND GENERALIZED TRAPPING STATES FOR STIRAP 

We now return to the case of ordinary STIRAP, i.e. without a detuning pulse D. The formal equivalence of W(í) 
and W(í) suggest an iteration of the procedure introduced in the last section. We define an nth order generalized 
adiabatic basis by the iteration: 



*< B >(t) 



= u r 



l(*) 4 



(n-l)(t) 
2 

(n-1), 



4> 



= U„_i(í)*·U„_ 2 (í)*···U (í) < 



1P1 

Í>3 



Correspondingly we obtain for the vector of probability amplitudes in the nth order basis 

BW = U„_iB("- 1 ) = U„_i • U n _ 2 ■ • ■ U C = V„C 
where we have dropped the time dependence. The nth order transformation matrix is defined as 



U n (í) 



1 

sin0„(í) 
i cos 0„ (í) 



with 



sin 0o (í) 
sin n (í) 



P{t) 

n (*)' 
n„_i(í) 



Q„(í) ' 

The iteration is illustrated in Figj|. 



cos 0o (í) 
cos 9 n (t) 



S(t) 
ílo(í)' 

20 n -l(í) 

n n (t) ' 





cos n (í) 

-i sin n (í) 



0„(í) = V"n-lW 2 +40n-l(í) 5 



(21) 

(22) 
(23) 

(24) 
(25) 
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FIG. 5. Iterative definition of nth order adiabatic basis 
In the roth-order basis, the equation of motion has then the form 



with 



W„(í) = - 



^BW(í) = -iW B (i)BW(í), 



Q n (i)sm0 n (í) 

Q n (i)sm0„(t) Í2 n (í)cos0„(t) 

n n (í)cos6»„(í) 

O n _i(í) 
íï»-i(*) 2è n _i(t) 
20„_i(t) 



(26) 



(27) 



If cos8k(t) vanishes, which implies that is time-independent, the state $3 decouples from the interaction. 

In this case exact analytic solutions of the atòmic dynamics can be found as discussed in the next section. The 
analytic solutions also include cases of population or coherence transfer. If cos 8k(t) does not vanish but is small, the 
corresponding coupling in the evolution matrix can be treated perturbatively. In such a situation we have a generalized 
adiabatic dynamics. 

In conclusion of this section it should be noted, that the iterative definition of a generalized adiabatic basis is 
conceptually very similar to the superadiabatic approach of Berry [ fi"H introduced for two-level systems. 



IV. GENERALIZED MATCHED PULSES AND ANALYTIC SOLUTION OF ATÒMIC DYNAMICS 



If a dynamical angle 9 n -i is a constant, the time-dependent state ^^\t) is decoupled from the interaction (constant 
of motion) . In this case the dynamical problem reduces to that of a two-state system interacting via a real resonant 
coherent coupling plus a decoupled state. 



d_ 
~dt 



B[ n \t) 



B, 



sy ; (í) 



o 

O n _i(í) 

o 



n n -x(t) o 
o o 
o o 



B[ n \t) 



This equation can immediately be solved 

B[ n) (t) = B[ n \0) cos^(í) - iB { 2 n) (0) sin^(í), 



BÍ n \t) = BÍ n) (0) cos^(í) 



iB[ n) (0) silici), 



4 B, (í) = C(o)> 



(28) 



(29) 
(30) 
(31) 



where 
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(32) 



In particular if the atom is initially in the trapping state, it will stay in that state. 

For example if 9o does not depenc! on time, the usual dark state is an exact constant of motion. As can be 
seen from Eq.(^), for 8q to be time- independent, Stokes and pump need to be either cw fields or need to have the 
same envelope function, i.e. have to be matched pulses [ p^[ , 

S(t) = O (í) cos O , (33) 
P(t) = fi„(i) sin0 O) (34) 

where íio(í) can be an arbitrary function of time and 6q =const. The atòmic dynamics is trivial in this case. Since 
is time-independent, the trapping state is a constant superposition of the bare atòmic states 1 and 3. 

On the other hand, if some higher-order dynamical angle 6 n is constant, the system remains in a generalized 
trapping state if initially prepared in it. The projection of this state onto the bare atòmic basis is in general time- 
dependent, and one can have a substantial rearrangement of atòmic level population including population transfer. If 
a higher-order dynamical angle is constant we will call pump and Stokes pulses generalized matched pulses. 

To obtain an explícit condition for generalized matched pulses in terms of P(t) and S(t) we successively integrate 
relations (|25|). This leads to the iteration 

Qk-i(t) = \ í dt'n k (t') co S e k (t') + 0° k , 

n k _ 1 (t) = n k (t)sme k (t), (35) 

starting with some 9 n {t) = n = const. and íín(í) as an arbitrary function of time. Each iteration leads to one constant 
k , which can be freely chosen. The application of generalized matched pulses to coherent population transfer will be 
discussed in the next section. 

As noted before there may be cases, where for some number n the dynamical angle 9 n (t) does depend on time 
but its time-derivative is much smaller than the corresponding generalized Rabi-frequency Q n (t), while the same is 
not true for all k < n. In this case the state $3 (t) is an approximate constant of motion and we have an nth 
order adiabatic process. The example of loop-STIRAP discussed in the last section is a realization of a higher-order 
adiabatic process, which is non-adiabatic in the first-order basis. 



V. APPLICATION OF GENERALIZED MATCHED PULSES TO POPULATION- AND COHERENCE 

TRANSFER 

In the following we discuss several examples for a coherent transfer of population from one non-decaying state to 
the other or to the excited state using generalized matched pulses. We furthermore discuss the possibility to transfer 
coherence, for example from the ground state transition to an optical transition. Since in all cases there exist a 
generalized trapping state which is an exact constant of motion, we can obtain exact analytic results for the atòmic 
dynamics. 



A. Population and coherence transfer with second-order generalized matched pulses 

1. Complete transfer of coherence from a ground-state doublet to an optical transition 
(2) 

First we discuss the case when $3 is an exact constant of motion, i.e. a trapping state. Furthermore we assume 
that the state vector coincides with this trapping state at t — —00. Then the system will remain in the trapping 
state at later times. Therefore B\ =const. and it is clear from Figj| that \& is a time independent superposition of 
states <&j and and thus has at all times a constant probability amplitude of the bare atòmic state 2. In fact 
from 



C = V^B^ 



v.; 



(36) 



8 



we find 



C 2 (í) 
Cz(t) 



-l COSt^i 



i tan Q\ cos 6*o (í) 
1 

—i tan 0i sin0o(t) 



(37) 



We now identify state 2 with a lower i.e. non-decaying level and state 3 with an excited states. The pump pulse 
P(t) then couples two ground states which could be realized for example by a magnètic coupling. The Stokes pulse, 
which couples states 2 and 3 is considered an optical pulse. Due to the finite and constant admixture of state 2 to 
the trapping state, second-order generalized matched pulses are best suited to transfer coherence for example from 
the 1-2 transition to the 3-2 transition. 

We now want to construct pulses, that would lead to the desired complete coherence transfer. To achieve this we 
have to satisfy the initial and final conditions 



o (-oo) = O, 
o (+oo) = tt/2. 

On the other hand, the iteration equation ([sH]) requires for second-order matched pulses that 

0o (*) - 5 / dt' fií (í') cos e 1 + e° , 

^ J -oc 

Q (t) = fií (f) sia 0i, 



(38) 
(39) 



(40) 
(41) 



where 0i and 0§ are arbitrary constants and fii(í) an arbitrary positive function of time. To fulfill the initial condition 
( |38| ) we set 0[J = 0. In order to satisfy the final condition (39) we then have to adjust the total pulse area (see Eq.(40)) 



/oo 
dífio(i) = 7rtan0 
-oo 

Thus pump and Stokes pulses have the form 

~ttA (t) 



P(t) = fi (í)sin 
S(t) = fi (í) cos 



2A . 
ttA (t) 



2A 



with 



A(t)= / dí'fio(í') 

J — oo 

With this choice an initial coherent superposition of states 2 and 1 

oo) = — i cos 0i ip2 + sin 0i ipi 
can be completely mapped into a coherent superposition of states 2 and 3 

^(+00) = —i cos 0i ip2 — sin 0i 1P3. 



(42) 

(43) 
(44) 

(45) 

(46) 
(47) 



In order to transfer a given ground-state coherence to an optical transition the pulse area Aq should be chosen 
according to J42]), A = \C\{— oo)/C2(— co)|. The shape of fi(í) is otherwise arbitrary. It should be noted that 
Eq.([46"|) requires a certain fixed phase of the initial coherent superposition. The phase of the pump pulse, which is 
included in the definition of tpi (cf. Sec. II), may need adjustment to satisfy this condition. 

In Fig.|| we have shown the populations of the bare atòmic states for the example fio(í) = \/tt exp(— i 2 ) (A = n) and 
^(—00) = l/v / 2('0i from a numerical solution of the Schròdinger equation. One clearly sees that all population 

from state 1 is transferred to state 3. This transfer happens without diabatic losses despite the fact that A = ir and 
thus the usual adiabaticity condition is only poorly fulfilled. 
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FIG. 6. Example of complete coherence transfer from 1 — 2 to 3 — 2 with second-order generalized matched pulses. Plotted 
are the populations of bare atòmic states for a pair of pulses as shown in the insert. One recognizes constant population in 
state 2 and complete transfer of population in 1 to 3. 

The process discussed here may have some interesting applications, since it allows to transfer coherence from a 
robust and long-lived ground state transition to an optically accessible transition. 

The population transfer from 1 to 3 with hnite constant state amplitude in 2 discussed here coincides with the 
solution found by Malinovsky and Tannor |]Ï3| with numerical optimization techniques. Assuming a hnite constant 
amplitude in state 2, these authors numerically optimized the peak Rabi-frequency (which in this case is the only 
remaining free parameter) to achieve maximum population transfer. They found, that in order to maximize the hnal 
amount of population in state 3, the peak Rabi-frequency has to be larger than a certain critical value. This can very 
easily be verihed from the generalized matched-pulse solutions ( pÏ2] , [47| ). 

|C3(oo)| 2 = S in 2 ^ = ^^ (48) 

|C 2 (oo)| 2 -cos 2 1 = ^- I . (49) 

In the limit 9\ — > tt/2, which implies A — > oo, the admixture of level 2 vanishes and we essentially have population 
transfer from state 1 to state 3. 




2. Population transfer from 1 to 3 and non-exponential diabatic losses 



We have seen in the last subsection that second-order matched pulses can be used to effectively transfer population 
from state 1 to 3, if there is an initial admixture of the excited state. This amplitude is inversely proportional to the 
square of the pulse area A. Therefore one could expect a good transfer for large A also if all population is initially 
in state 1. In this case there is some hnite amount of population which is not trapped in the generalized dark state 
$à . Clearly in order to achieve maximum population transfer, pump and Stokes pulse should be in countcrintuitive 

Since the pulses are assumed to be second-order matched 



order and hence conditions (|38|) and ( |39| ) should be fulhlled. 
pulses, the dynamical problem with the initial condition 



\B{ 2 ' 


(-oc) 




"1" 









(-oo) 


= Ui-U 







i cos 9i 


B?' 


(-oo) 









sin 6-y 



(50) 



can easily be solved (see Eq.(|^]3|)). From Eqs.@ we find ü x {t) = Q (*)/ sin 6 X . Thus 

\^TT 2 +A 2 , 



Bi 2) (oo) 



B«(oo)=t 



Vtt 2 + A 2 



sm 



fíf (oo) 



Vtt 2 + A 2 
A 



cos 



-\Ztt 2 + A 2 



Vtt 2 + A 2 ■ 



(51) 
(52) 
(53) 
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where A is the total pulse area defined in (|42|). From this we find the asymptotic populations of the bare atòmic 
states 



Ci(oo) 



c 2 (oo; 



C 3 (oo) 



— pr sin 2 (\ \/n 2 + A 2 



A 2 
4ir 2 A 2 



(tt 2 + A 2 ) 
1 



sin 4 (1^ + aA , 
A 2 + 7r 2 cos | i Vti -2 + -4 2 



-i 2 



(tt 2 + A 2 ) 

Thus the diabatic losses scale in general with 1/A 2 , i.e. non-exponentially with A. Furthermore for 

l;Vn 2 + A 2 = 2n-K or A = n\J 16ti 2 - 1 



(54) 
(55) 

(56) 
(57) 



with n = 1,2,... the population transfer is complete. We show in Fig. |^ the final population in state 3 as a function 
of A/ir. 
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15 



20 



FIG. 7. Final population in state 3 as a function of total pulse area A/n for population transfer from state 1 with second-order 
matched pulses. For A/n — \J16n? — 1 the transfer is complete (100.00%) 

A special case of the population transfer with second-order matched pulses discussed in the present section is the 
analytical model discussed by Vitanov and Stenholm in These authors considered a pulse sequence with 

íïo(í) = ^sech 2 (Çj , O (t) = j tanh(~)+l (58) 

and thus tan é'l = Qo(í)/26'o(t) = a/n = const. 

B. Population transfer via large-area third-order matched pulses 

Next we analyze the possibility of population transfer when $3 is exactly trapped. In order for $ 3 to be a 
constant of motion or equivalently to have third-order matched pulses #2 =const. We assume again that the system 
state vector ^ is initially in the trapping state in which it will remain for all times. In order to realise population 
transfer from state 1 to state 2 or 3 in this case, we furthermore must satisfy the initial conditions 



d (-00) = 1, C 2 (-00) = 0, C 3 (-00) = 0. 



(59) 



This can be translated into a condition for the initial vàlues of the dynamical phases #0 and 6\ using Eq.(g2|). In fact 
from 



V3 W = V3- 1 



(60) 
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we find 



Ci (i) = i (cos O (í) sin 0i (í) sin 2 - sin O (í ) cos 2 ) (61) 
C 2 (í) = -i cos 0i (í) sin 2 (62) 
C 3 (t) = -i (cos0q(í)cos0 2 +sm0 o (í)sin0i(í)sin0 2 ) . (63) 



The initial condition it is fulfillcd whcn 



cos 0o (— oo ) sin 0i (— oo) sin 02 — cos0 2 sin0 o (— oo) =1, (64) 

sin 2 cos 0i (— oo) =0, (65) 
cos 0o (—oo) cos 02 + sin 02 sin 0q (—oo) sin 0i (—oo) =0. (66) 



The result is 

7T „ „ . 7T 



1 (-oo) = ~, 0o(-oo) = 2 + -. (67) 
Erom Eq.(poj) we find the following differential equation 

2 d0 ^ = afii(í), where a = (tan0 2 )" 1 = const. (68) 

Introducing 

x(t) =tan0i(í) (69) 

we find furthermore 

2x x 



{í + x 2 ) 3/2 

de (t) _ a (t) 



dt x{t) 

Integrating these equations and taking into account the initial conditions (J67j) yields 



aíl (t), (70) 

(71) 



lan0,(/l= /■(/) = 1- ''- , ) 



7T 1 



where 



e (t) = e 2 + - + - 1-^/1-/2(4) , (73) 

2 a L J 



/(*) = ? /' dí'^oCí')- (74) 



2 



— oo 



f?o(í) is an arbitrary smooth function which we assume to vanish at infinity, íío (i 00 ) — 0. We still have one free 
constant a, which we can choose. As we will show now, we can choose a such that the efficiency of the transfer from 
statc 1 to states 3 or 2 approaches unity. 

1. Population transfer from ground state to state 3 

In order to transfer the initial population from state 1 to the target state 3, it is necessary to satisfy the final 
conditions 

0i(+°°) = |, 0o(+^)=^2 (75) 

which implies 
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tan 0i (oo) 
o (oc) 



4 






- 1 


a 2 A 2 


7T 


1 " 


+ 2 4 


o . 



oo, 



where 



A 



1 - ^/l-a 2 A 2 /A 



dtü (t) 



is the pulse area. From these condition one finds the constraint 

47T 



A 2 - 



The diabatic losscs in the limit A ^> 1 are 



1 



C 3 (oo) 



A > 7T. 



4^ 2 
2tt 2 + A 2 



(76) 
(77) 

(78) 
(79) 
(80) 



and thus in the adiabatic limit we have essentially complete population transfer from state 1 to state 3. 

Fig.|| shows an example of population transfer with third-order matched pulses. Here íïo(t) = ^4/2sech 2 (í) and 
A = 20n. Pump and Stokes pulses are shown in the upper frame and the population histories in the lower one. We 
see that the amplitudes of the Stokes and pump pulses are unequal. As in ordinary STIRAP the population of the 
state 2 is small during the evolution. 





2(x 10) 



-2-101234 
t 

FIG. 8. Population transfer from 1 to 3 with third-order matched pulses. Upper frame shows pulses, lower frame population 
dynamics. Here fio (í) = j4/2sech 2 (í) and A/ty — 20. 



2. Population transfer from ground state to the state 2 
In order to transfer the initial population from state 1 to state 2, it is necessary to satisfy the conditions 

e x (+oo) = o, e 2 = | (8i) 



In this case we have to fix a to be 



a=j, A»l. (82) 
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Fig.^ shows the pulses P(t) and S(t) and the evolution of the atòmic populations. Here fio(í) = ^4/2sech 2 (í) and 
A = 16. We see that the Stokes and pump pulses are in a counterintuitive sequence. At first the atòmic population 
oscillates between state 1 and 3, but as the pulse sequence proceeds the whole population is transfered into state 2. In 
other words, during the full pulse sequence there occur several STIRAP transitions, but due to the large nonadiabatic 
coupling the population accumulatcs in state 2. 




-2-10 1 2 3 4 
t 

FIG. 9. Population transfer from 1 to 2 with third-order matched pulses. Upper frame shows pulses, lower frame population 
dynamics. Here ílo{t) = A/2 sech 2 (í) and A = 16. 



VI. SUMMARY 



We have introduced the concept of generalized dressed states in order to explain the success of population transfer 
in stimulated Raman adiabatic passage with a loop coupling. If the interaction of a three-level system with a pair of 
time-dependent pump and Stokes pulses is described in terms of the so-called dark and bright states instead of the 
instantaneous eigenstates of the Hamiltonian, the original three-state-two-ficld system is transformed into a system 
of three states coupled by two effective interactions ^jQ- This allows for an iteration procedure leading to higher- 
order adiabatic basis sets O]. We showed that in the case of loop-STIRAP there is a higher-order trapping state, 
which is an approximate constant of motion even when the usual adiabaticity condition is not fulfilled. This state 
adiabatically rotates from the initial to the target quantum state of the atom and thus leads to efficient population 
transfer, however, at the expense of placing some population into the decaying atòmic state. 

The concept of generalized trapping states allows the construction of pulse sequences which lead to an optimum 
population or coherence transfer also for small pulse areas and allows for solutions for the atòmic dynamics. If pump 
and Stokes pulses fulfill certain conditions (so-called generalized matched pulses) the effective 3x3 coupling matrix 
factorizes at a specific point of the iteration. The trapping state of the corresponding nth-order adiabatic basis is then 
an exact constant of motion. In this case the atòmic dynamics reduces to a two-level problem with a real coupling 
which can be solved analytically. 

For ordinary matched pulses, i.e. if pump and Stokes have the same shape, the atòmic dynamics is rather limited. 
The corresponding dark state is a constant superposition of states 1 and 3. In the case of generalized matched pulses, 
however, the trapping state has a time-dependent overlap with the bare atòmic states and thus population or coherence 
transfer is possible. We have discussed with specific example population transfer with second and third-order matched 
pulses. We found that for certain vàlues of the pulse areas complete population or coherence transfer is possible. In 
the general case the diabatic losses scale non-exponentially with the inverse pulse area. 
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